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We present a Dalitz plot analysis of chaxmless decays to the final state tt^tt^tt^ using a 
sample of (465 ± 5) x 10® BB pairs collected by the BABm experiment at -v/i = 10.58 GeV. We 
measure the branching fractions B{B^ tt^-k^tt^) = (15.2 ± 0.6 ± 1.2 ± 0.4) x 10"'', B{B^ 
p°(770)7r±) = (8.1 ±0.7± 1.21?:*) x 10"^ B(B^ /2(1270)7r±) = (1.57±0.42±0.16 -^111) x 10-^ 
and B{B^ tt^tt^tt^ nonresonant) = (5.3 ± 0.7 ± 0.61q 5) x 10~®, where the uncertainties are sta- 
tistical, systematic, and model-dependent, respectively. Measurements of branching fractions for the 
modes B^ — » p''(1450)7r^ and B^ — > /o(1370)7r* are also presented. We observe no significant di- 
rect CP asymmetries for the above modes, and there is no evidence for the decays B^ — » /o(980)7r*, 
B"^ XcoTT*, or B^ Xc27r*. 
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INTRODUCTION 

Decays of B mesons to three-body charmless fi- 
nal states probe the properties of the weak interac- 
tion through their dependence on the complex quark 
couplings described in the Cabibbo-Kobayashi-Maskawa 
(CKM) matrix P, Q ■ Furthermore, these decays test dy- 
namical models for hadronic B decays. 

One can measure direct CP asymmetries and constrain 
magnitudes and phases of the CKM matrix elements 
using individual channels that appear as intermediate 
resonances in the B^ — > tt^tt^tt^ decay. For example, 
the CKM angle 7 could be extracted from the interfer- 
ence between the decay B^ Xco^^^, which has no CP- 
violating phase (in the standard paramctrization), and 
other modes such as B^ p"(770)7r± [3, 4, 5, gM^- 



Studies of B^ 



r±7r±7rT 



can also be useful for a pre- 



cise measurement of the CKM angle a. A theoretically 
clean determination of this angle can be obtained from 
the decay-time dependence of the interference between 
B° ^ p+TT", B° ^ /9"7r+, and B° p°7r° via the anal- 
ysis of the Dalitz plot for B^ — > 7r+7r~7r° deca ys i9|| (r e- 
cently implemented by BaBAR [13] and Belle [liriil)- 
Charged B decays offer a large statistics sample with 
which to determine additional resonant or nonresonant 
contributions to the three-pion Dalitz plot that can af- 
fect the measurement of a. For example, the Dalitz plot 
analysis of B^ tt^tt^tt^ allows one to check for effects 
from B^ Lj{782)'K^, that could cause large direct CP 
violation due to p-uj mixing 13]. It is particularly impor- 
tant to limit the possible effects of broad scalar structures 
[including the so-called /o(600) or a] and nonresonant 
contributions [H [H 0, [13, [11, [11] . 



Furthermore, a number of unexplained structures have 
been observed in charmless B decays to Ktttt 2C , 21 , 

'MM, m\ 



22, 23, 24, \M, KKtt m, \M, and KKK 



final states. Verifying the presence of these structures 
m — > TT^TT^TT^ decays would help to determine their 
nature and involvement in hadronic B decays. 

In this paper we present an amplitude analysis of 
B^ TT^n^n^ decays based on a 424 fb~^ data sample 
containing (465 ± 5) x 10^ BB pairs iN_^). The data 
were collected with the BaBAR detector [3C ]at the PEP- 
II asymmetric-energy e^e^ storage rings 3l| operating 
at the T{AS) resonance with center-of-mass (CM) energy 
of y/s = 10.58 GeV. An additional total integrated lumi- 
nosity of 44 fb"^ was recorded 40MeV below the T{AS) 
resonance ( "off-peak" data) and was used to study back- 
grounds. Compared to our previous publication [s^], in 
addition to doubling the data sample we have included 
several improvements in reconstruction algorithms that 
enhance the signal efficiency, made numerous modifica- 
tions to the analysis to increase the sensitivity to direct 
CP violation effects (for example, by including more dis- 
criminating variables in the maximum likelihood fit), and 



improved our model of the Dalitz plot structure. 

The remainder of the paper is organized as follows: 
Sec. II describes the amplitude analysis formalism. Sees. 
Ill and IV give details about the selection of signal B 
decays and how backgrounds are considered. Sec. V 
presents the results from the likelihood fit. Sec. VI gives 
an account of the various sources of systematic uncer- 
tainties, while Sec. VII summarizes the results. 



AMPLITUDE ANALYSIS FORMALISM 



A number of intermediate states contribute to the de- 
cay B^ Tr^TT^TT^. We determine their contributions 
with a maximum likelihood fit to the distribution of 
events in the Dalitz plot. This procedure has been de- 
scribed in detail in our previous publications [23,[2l|,[3l. 

The B^ TT^TT^TT^ decay contains two same-sign pi- 
ous in the final state. We distinguish these particles ac- 
cording to the invariant mass they make when combined 
with the oppositely charged pion, and draw the Dalitz 
plot in terms of heavy and light invariant masses-squared 
of the tt'^tt^ systems (denoted m^^x ^^'^ "^min' respec- 
tively), so that each candidate has a uniquely defined po- 
sition. Moreover, we explicitly enforce the symmetriza- 
tion of the total amplitude under exchange of identical 
bosons. 

The total signal amplitudes for B^ and B^ decays are 
given by 

^ = ^(^max; '^min) = / y "^j^j ("^max' ^inin) ' (1) 



^ = ^(™max>"lmm) J^'^^^- 



' J (^'^maxJ '^min) 



The complex coefficients Cj and Cj for a given decay mode 
J contain all the weak phase dependence. Since the Fj 
terms contain only strong dynamics, Fj = Fj. We use 
the following paramctrization [2]| for the amplitude co- 
efficients: 



(.Xj + Aa;^) +i{yj + Ay^) 
{xj - Axj) + i{yj - Ayj) . 



(2) 



In this approach, Xj and yj [Axj and Ayj) are the CP- 
conserving (-violating) components of the decay ampli- 
tude. 

The Fj distributions describe the dynamics of the de- 
cay amplitudes and are written as the product of an in- 
variant mass term Rj , two Blatt-Weisskopf barrier form 
factors Xj, and an angular function Tj 

^j(n^max, "imin) = R'j{rn) ^jip*) ^jil) Tj{m) , (3) 

where m (J) is the mass (spin) of the resonance, p* is the 
momentum the bachelor pion that is not part of the reso- 
nance in the B meson rest frame, and q is the momentum 
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of either daughter in the rest frame of the resonance (we 
use the c = 1 convention for all equations in this paper). 
The Fj are normalized over the entire Dalitz plot: 



\Fj{m 



2 

max 



1 '""min 



) I ^ dm 



2 

max 



(4) 



The Blatt-Weisskopf barrier form factors [33| are given 



by: 



Xj^oiz) = 1, 

Xj=i{z) = Vl/il + izTBwV], 

Xj^2{z) = v/l/[(^^Bw)^ + 3(zrBw)' + 9]^ 



(5) 



where the meson radius parameter tbw is taken to be 
4.0 ± 1.0 (GeV/c)-i [33. 

For most resonances in this analysis the Rj are taken 
to be relativistic Breit-Wigner lineshapes 

-RjM = 2V~- "FT^ ' 

(rriQ — TO^j — imoL (m) 

where mo is the nominal mass of the resonance and T{m) 
is the mass-dependent width. In the general case of a 
spin- J resonance, the latter can be expressed as 



/ m Xjiqo) ' 

The symbol Tq denotes the nominal width of the reso- 
nance. The values of mg and Tq are obtained from stan- 
dard tables 3J] when they are well known. The symbol 
qo denotes the value of q when m ~ toq- 

The angular distribution terms 7} in Eq. ([3]) follow 



35l . l36l | . For the decay of 



the Zemach tensor formalism 
a spin zero B-meson into a spin J resonance and a spin 
zero bachelor particle this gives [s^l 



,7=0 



.7=1 



1, 



(8) 



~2p-q, 

\[3ip-q)^-i\p\\q'\f 



where p is the momentum of the bachelor particle and q 
is the momentum of the resonance daughter with charge 
opposite from that of the bachelor particle, both mea- 
sured in the rest frame of the resonance. 

The Gounaris-Sakurai parametrization (ssj of the P- 
wave scattering amplitude for a broad resonance decay- 
ing to two pions is used for the p^{770) and p°(1450) 
lineshapes 

1 4- To d/mo 



Rj (to) 



(toq — TO^) -I- /(to) — i TOor(m) 



where 
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and the function h{m) is defined as 
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The normalization condition at Rj (0) fixes the parameter 
d — /(0)/(roTOo). It is found to be 



7^ 9o 



In 



2m,r 



^0 



(13) 



We model the nonresonant component using an empir- 
ical function that has been found to accurately describe 
nonresonant contributions in other charmless three-body 
B decays [2i[24, 25, 28, 29]: 



^). 



(14) 



We include this term in the coherent sum given by Eq. ([T]) 
when calculating the total signal amplitude over the 
Dalitz plot. 

To allow comparison among experiments we present 
results also in terms of fit fractions {FFj), defined as the 
integral of a single decay amplitude squared divided by 
the total matrix element squared for the complete Dalitz 
plot 
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(15) 



Note that the sum of all the fit fractions is not neces- 
sarily unity due to the possible presence of constructive 
or destructive interference. The CP asymmetry for each 
contributing resonance is determined from the fitted pa- 
rameters 



A 



CP, J 



(16) 



~2{xjAxj +yjAyj 



(x,)2 + {Axjf + {yj)^ + {Ay,y 



The signal Dalitz plot probability density function 
(PDF) is formed from the total amplitude as follows: 



Psig{m 
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where qs is the charge of the i?-meson candidate, and 



e{m 



max' ' mm 



) and e = e{m 



max' "mm 



) are the signal 



reconstruction efficiencies for and B events, respec- 
tively, defined for all points in the Dalitz plot. 
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CANDIDATE SELECTION 

We reconstruct B candidates from events that have 
four or more charged tracks. Each track is required to 
be well measured and to originate from the beam spot. 
They must have a minimum transverse momentum of 
50 MeV/c, and a distance of closest approach to the beam 
spot of less than 1.5 cm in the transverse plane and less 
than 2.5 cm along the detector axis. B candidates are 
formed from combinations of three charged tracks, and 
particle identification (PID) criteria are applied to reject 
electrons and to separate pions from kaons. In our fi- 
nal state, the average selection efficiency for pions that 
have passed the tracking and PID requirements is about 
93% including geometrical acceptance, while the average 
misidentification probability of kaons as pions is close to 
8%. 

Two kinematic variables are used to identify signal B 
decays. The first variable is 



(18) 



the difference between the reconstructed CM energy of 
the S-meson candidate (E^) and \/s/2, where ^/s is 
the total CM energy. The second is the beam-energy- 
substituted mass 



^Es = y s/A- 



(19) 



where p*^ is the B momentum measured in the CM frame. 
The Toes distribution for signal events peaks near the 
B mass with a resolution of around 2.5MeV/c^, while 
the AE distribution peaks at zero with a resolution of 
approximately 20MeV. We initially require events to 
lie in the region formed by the following selection cri- 
teria: 5.200 < TOES < 5.286 GeV/c2 and -0.075 < AE < 
0.300 GeV. The region of AE below -0.075 GeV is heav- 
ily contaminated by four-body B decay backgrounds and 
is not useful for studying the continuum background. 
The selected region is then subdivided into three ar- 
eas: the "left sideband" (5.20 < toes < 5.26GeV/c2 
and |Ai5| < 0.075 GeV) used to study the background 
AE and Dalitz plot distributions; the "upper sideband" 
(5.230 < TOES < 5.286 GeV/c^ and 0.1 < AE < 0.3 GeV) 
used to study the background toes distributions; and 
the "signal region" (5.272 < toes < 5.286 GeV/c^ and 
|A£:| < 0.075 GeV) with which the final fit to data is 
performed. Following the calculation of these kinematic 
variables, each of the B candidates is refitted with its 
mass constrained to the world-average value of the B 



meson mass |34| in order to improve the Dalitz plot posi- 
tion resolution and to make sure all events lie within the 
kinematic boundary of the Dalitz plot. 

The dominant source of background comes from light- 
quark and charm continuum production (e+e^ qq, 
where q = u,d,s,c). This background is suppressed 
by requirements on event-shape variables calculated in 



the CM frame. We compute a neural network (NN) 
from the following five variables: the ratio of the second- 
and zeroth-order angular moments {L2/Lq), with Lj = 
^ j Pi I cos 9i p , where 9i is the angle of the track or neutral 
cluster i with respect to the signal B thrust axis, pi is its 
momentum, and the sum excludes the daughters of the 
B candidate; the absolute value of the cosine of the angle 
between the direction of the B and the detector axis; the 
magnitude of the cosine of the angle between the signal 
B thrust axis and the detector axis; the output of a mul- 
tivariate i3-fiavor tagging algorithm |39| multiplied by 
the charge of the B candidate; and the ratio of the mea- 
sured proper time difference of the two B decay vertices 
and its statistical uncertainty. We train the NN using 
samples of off-peak data and signal Monte Carlo (MC) 
events generated with the phase-space distribution. A 
selection requirement is imposed on the NN output that 
accepts about 48% of signal events while rejecting 97% 
of continuum background events. 

Dalitz plot distributions of the reconstruction effi- 
ciency for and B^ events are modeled with two- 
dimensional histograms formed from a sample of around 
7 X 10'^ B"^ -K^TT^-K^ phase-space MC events. All se- 
lection criteria are applied except for the exclusion of 
certain invariant-mass regions described below. We take 
the ratio of two histograms, the denominator contain- 
ing the true Dalitz plot distribution of all generated MC 
events and the numerator containing the reconstructed 
MC events. The reconstructed events are weighted in or- 
der to correct for differences between data and MC sim- 
ulations in the tracking and PID efficiencies. In order to 
give better resolution near the edges of the Dahtz plot, 
where most reconstructed events lie, the histograms are 
formed in the "square Dalitz plot" [13, [sil coordinates. 
We use 50x50 bins and smooth these histograms by 
applying linear interpolation between neighboring bins. 
The efficiency shows very little variation across most of 
the Dalitz plot but decreases towards the corners where 
one of the particles has low momentum. The effect of 
experimental resolution on the signal model is neglected 
since the resonances under consideration are sufficiently 
broad. The average reconstruction efficiency for events in 
the signal region for the phase-space MC sample is about 
15%. The fraction of misreconstructed signal events is 
only 5%, and MC studies indicate that there is no need 
for any explicit treatment of these events. 



BACKGROUNDS 

In addition to the continuum {qq) background wc also 
have backgrounds from BB events. There are four main 
sources: (i) combinatorial background from three un- 
related tracks; (ii) three- and four-body B decays in- 
volving an intermediate D meson; (iii) charmless two- 
and four-body decays with an extra or missing parti- 
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cle; and (iv) three-body decays with one or more parti- 
cles misidentified. We reject background from two-body 
decays of D mesons and charmonium states by exclud- 
ing invariant masses (in units of GeV/c^) in the ranges: 



1.660 < 



< 1.920, 3.051 < 



< 3.222, 



and 3.660 < to^+tt- < 3.820. These ranges reject de- 
cays from D'^ ^ K^Tr~ (or tt+tt^), J/ip , and 
■0(25) t^t~ respectively, where ^ is a lepton that 
has been misidentified as a pion. We also employ a 
special requirement to reject the decay process — *■ 
K^gV^\K^^ 7r+7r~, by excluding candidates where the 
vertexed mass of two oppositely charged pions lies in the 
range of [478, 516] MeV/c^. 

We use a large sample of MC-simulated BB decays, 
equivalent to approximately 3 times the integrated lu- 
minosity of the data sample, to identify the important 
B backgrounds that survive the invariant-mass exclusion 
requirements described above. In total, 53 i3-meson de- 
cay modes are identified for which larger samples of ex- 
clusive MC events are used for further study. We com- 
bine modes that have similar behavior in the discrimi- 
nating variables tties and into a B-background cat- 
egory. There are four such categories: the first contains 
the two-body decays B^ tt+tt^ and 5° K~tt^ , the 
second is dominated by B^ —> K^tt^tt^ and contains 
other decays with similar topologies, the third contains 
only — > 7r+7r~7r'^, and the fourth contains the remain- 
ing backgrounds from B decays that are combinatorial 
in nature. For each i?-background category the com- 
bined mESj A_E, and Dalitz plot distributions are created 
where the relative contributions of various decay modes 
in a specific category are calculated from the reconstruc- 
tion efficiencies from MC simulations and the branching 
fractions listed by the Particle Data Group |3^] and the 
Heavy Flavor Averaging Group [i^] ■ These distributions 
are used in the likelihood fit described below. 

Background Dalitz plot distributions arc included in 
the likelihood fit through the use of two-dimensional 
histograms. For backgrounds from B decays these his- 
tograms are formed from the various MC samples. For 
the continuum background the left sideband data sam- 
ple is used. Since this data sideband also contains events 
from B decays, MC samples are used to subtract these 
events. To these B-subtracted sideband events, we add 
off-peak data events from across the whole range of toes 
and AE in order to enhance statistics. We have ver- 
ified that the shapes of various discriminating variables 
are compatible between the sideband and off-peak events. 
As for the reconstruction efficiency histograms, the back- 
ground Dalitz plot distributions are formed in the square 
Dalitz plot coordinates and are smoothed by linear in- 
terpolation applied between neighboring bins. Separate 
histograms are constructed for B~^ and B~ events. The 
qq- and i?-background PDFs are identical in their con- 



struction, and the qq PDF is shown here as an example: 
1 
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where Aqq is the charge asymmetry in the background, 
and (5(TO^a^, TO^i„) and Q{m'^,a^^,m^^^) are the Dahtz 
plot distributions of qq events in selected B^ and B^ 
samples, respectively. 



MAXIMUM LIKELIHOOD FIT 

To provide further discrimination between the signal 
and background hypotheses in the likelihood fit, we in- 
clude PDFs for the kinematic variables toes and Ai5, 
which multiply that of the Dalitz plot. The signal toes 
shape is modeled with the sum of a Gaussian function 
and a Crystal-Ball lineshape 



41| 



and the AE shape is 
modeled with a double Gaussian function. The param- 
eters of these functions are obtained from a sample of 
B^ Tr^TT^Tr'f MC events, modeled according to the 
Dalitz plot distribution from Ref . , and are appropri- 
ately adjusted to account for possible differences between 
data and MC simulations determined with a control sam- 
ple of £?+ D^Ti^-. K^TT^ decays. These param- 
eters are fixed in the fit to data. 

The qq toes distribution is modeled with the experi- 
mentally motivated ARGUS function [i^l . The end point 
for the ARGUS function is fixed to 5.289 GeV/c^, and the 
parameter describing the shape is fixed to the value de- 
termined from the combined sample of upper sideband 
and off-peak data. We model the continuum AE shape 
using a linear function, the slope of which is fixed to 
the value determined from the left sideband and off-peak 
data. The BB background distributions are modeled 
with histograms obtained from the mixture of BB MC 
samples. The yields of signal and qq events are allowed 
to vary in the final fit to the data while the yields of 
BB backgrounds are fixed to 11 (two-body decays), 195 
(B± ^ K^TT^TT^ type), 117 (5° ^ vr+Tr-vr"), and 495 
(combinatorial) events. 

The complete likelihood function is given by: 



(21) 
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NkViiml^^^, TO^i„, TOES, Ai;, qs) 
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where N is equal to Nk, Nk is the yield for the event 
category k, is the total number of events in the data 
sample, and 7^^ is the PDF for the category A; for event 
J, which consists of a product of the Dalitz plot, tties, 
and Ai? PDFs. The function — 21n>C is minimized in an 
unbinned fit to the data. 

Our nominal signal Dalitz plot model comprises a 
momentum-dependent nonresonant component and four 
intermediate resonance states: yO°(770)7r='=, /9'^(1450)7r='=, 
/2(1270)7r^, and /o(1370)7r^. The parameters used to 
describe these states are summarized in Table HI We fit 
4335 B candidates in the signal region selected from the 
data to obtain the central values of the Xj , Axj , j/j , and 
Aj/j parameters for each component, and use Eqs. pS]) 
and (Uni) to calculate the fit fractions and CP asymme- 
tries. We use p*'(770)7r='= as the reference amplitude, fix- 
ing its X, y, and Ay parameters to unity, zero, and zero, 
respectively. The signal yield, qq background yield and 
asymmetry are also free parameters of the fit, giving a 
total of 20 free parameters. 

The Dalitz plot model was determined using the re- 
sults of our previous analysis [32*1 and the changes in the 
fit likelihood and values when omitting or adding res- 
onances. The latter is calculated from the projection of 
the fit results onto the Dalitz plot using the formula 



X 



E 

1=1 



[y^ - fi^i)? 



(22) 



where yi is the number of data events in bin i and f{xi) 
is the number of events in that bin as predicted by the fit 
result. The number of degrees of freedom is calculated as 
rUf — h—l, where rib is the total number of bins used and h 
is the number of free parameters in the fit. A minimum of 
20 entries in each bin is required; if this requirement is not 
met then the neighboring bins are combined. Typically, 
Ub takes values around 100. 

In our previous study we found significant contribu- 
tions from pO(770)7r± and /2(1270)7r±; with /o(980)7r±, 
p'^(1450)7r^, and a uniform nonresonant term also in- 
cluded in the model. Because of the larger data sam- 
ple and many improvements to the analysis, we find 
it necessary to include an additional contribution from 
/o(1370)7r^, and to use a momentum-dependent nonres- 
onant amplitude [see Eq. (|14p ] in order to achieve a rea- 
sonable agreement of the fit with the data. We do not 
find any significant signal from /q (980)7r* , so we exclude 
this channel from our nominal model and calculate an up- 
per limit for its fit fraction. The statistical significance of 
the presence of a component is estimated by evaluating 
the difference A In £ between the negative log- likelihood 
of the nominal fit and that of a fit where all of the x, 
y, Ax, and Ay parameters for the given component are 
fixed to zero. This is then used to evaluate a p value 



P 



f{z;nd) dz , 



(23) 



where f{z; Ud) is the PDF of the distribution and Ud 
is the number of degrees of freedom, four in this case. 
We then determine the equivalent one-dimensional sig- 
nificance from this p value. We find that the /2(1270) 
contribution has a statistical significance of G.lcr, the 
p0(1450) 4.6a and the /o(1370) 3.9cr. 

Since the mass and width of the /o(1370) state are not 
well known [s^ ]. we determine the preferred values from 
data by scanning the likelihood values obtained with dif- 
ferent parameters. The mass and width are determined 
to be m^,(i37o) = 1400 ± 40MeV/c2 and r;„(i37o) = 
300 ± 80MeV, with a correlation of (-39 ± 4)%, where 
the errors are statistical only, and are obtained from a fit 
to the two-dimensional likelihood profile. Similarly, we 
determine the parameter of the nonresonant lineshape to 
be Qfnr — 0.28 ± 0.06 GeV~^c'' (statistical uncertainties 
only). 

Possible contributions from Xco""^ and Xc2T^^ are not 
significant so we set upper limits on their branching frac- 
tions. Furthermore, we do not find any evidence for a 
very broad enhancement at low tt^tt^ invariant mass 
such as could be caused by the decay an^. 

TABLE I: Parameters used to describe intermediate states 
in our nominal model. GS and RBW refer to the Gounaris- 
Sakurai and relativistic Breit-Wigner lineshapes, respectively. 

Resonance Lineshape Mass (MeV/c^) Width (MeV) 
p"(770) 



Ref. 



p°(1450) 
/2(1270) 
/o(1370) 



GS 
GS 
RBW 
RBW 



775.49 ± 0.34 

1465 ± 25 
1275.1 ± 1.2 
1400 ± 40 



149.4 ±1.0 [34] 
400 ±60 [34| 

185.0^21 [34] 

300 ± 80 See text 
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FIG. 1: (Top) signal and (bottom) qq distributions of (left) 
mES and (right) AE obtained from the fit to data using event- 
by-event signal and qq background probabilities 14i||. The 
solid lines show the PDF shapes used in the fit. 



Figure [T] shows the toes and AE distributions of signal 
and qq background determined from the fit with event- 
by-event signal and qq background probabilities for each 
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FIG. 2: Background-subtracted Dalitz plot of the combined 
— > n^Tv^n^ data sample in the signal region. The plot 
shows bins with greater than zero entries. The area of the 
boxes is proportional to the number of entries. The depleted 
bands are the charm and charmonia exclusion regions. 



candidate event [43| . The background-subtracted Dalitz 



plot of the data in the signal region can be seen in Fig. [2l 
The per number of degrees of freedom of the projec- 
tion of the fit result onto the Dalitz plot is 82/84. Using 
the fitted signal distribution we calculate the average re- 
construction efficiency for our signal sample to be 18%. 

We generate a large number of MC experiments with 
the fitted parameters, and from the spread of results of 
fits to those experiments we determine the statistical un- 
certainties on the parameters, FFj, and Acpj- This 
procedure takes into account correlations between the 
Xj, Axj, Uj, and Ayj parameters. The linear correla- 
tion coefficients between the FFj and Acp, j parameters 
are also obtained and are presented in Appendix . In 
order to calculate the branching fraction for an interme- 
diate mode, we multiply the fit fraction of the latter by 
the total inclusive tt^tt^tt^ branching fraction. 

They are needed for comparison with previous measure- 
ments and theoretical predictions. The tt^tt^tt^ signal 
yield is found to be 1219 ± 50 ± 75^24 events and the 
inclusive CP asymmetry to be (-1-3.2 ± 4.4 ± 3.1 %, 
where the uncertainties are statistical, systematic, and 
model-dependent, respectively. Additionally, the total 
yield and CP asymmetry of the continuum background 
are found to be 2337 ± 62 events and {+0.2 ± 2.7)%, 
respectively, where the uncertainties are statistical only. 
Further results are shown in Tables |TT] and IIIII 

Projections of the data, with the fit result overlaid, as 
TT^TT^ invariant-mass distributions can be seen in Fig. [3l 
A detailed examination of possible direct CP violation 
effects in the low tt^tt^ invariant-mass region is shown 
in Fig.m where we have subdivided the data into positive 
and negative values of cos 9h — P'Q /{IpWqI), where Oh is 
the helicity angle, pis the momentum of the bachelor par- 



ticle and q is the momentum of the resonance daughter 
with charge opposite from that of the bachelor particle, 
both measured in the rest frame of the resonance. The 
agreement between the fit result and the data is generally 
good; the per number of non-zero bins for these plots 
varies between 35/46 and 34/24. 

We calculate 90% confidence-level (CL) upper limits 
for components not included in the nominal Dalitz plot 
model. These are obtained by generating many MC ex- 
periments from the results of fits to the data where the 
extra component is added to the nominal Dalitz plot 
model, with all major systematic sources varied within 
their la uncertainties. We fit these MC samples and 
plot the fit fraction distributions. The 90% CL upper 
limit for each fit fraction is the value which includes 90% 
of the MC experiments. The branching fraction upper 
limit is then the product of the fit fraction upper limit 
and the total branching fraction for tt^tt^tt^ . 
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FIG. 3: (color online) Dipion invariant mass projections: 
(left) in the p°(770) region; and (right) in the regions of XcO 
and Xc2- The data are the points with statistical error bars, 
the dark-shaded (red) histogram is the qq component, the 
light-shaded (green) histogram is the BB background contri- 
bution, while the upper (blue) histogram shows the total fit 
result. The dip near 0.5 GeV/c^ in the left plot is due to the 
rejection of events containing Ks candidates. 



We have searched for the presence of multiple solu- 
tions in the fit to data with the nominal model. We find 
a second solution with a value of —2 In £ about ten units 
higher than our nominal fit, and with a of the Dalitz 
plot projection increased by four units. A comparison 
of the results between the two solutions is given in Ap- 
pendix . The most significant difference is seen in the 
/o(1370) fit fraction, which is much smaller in the second 
solution. In Dalitz plot analyses of Ktttt and KKK final 
states [13, [m, m, m [iH, m, nil , slmUar phenomena 
relating to multiple solutions have been observed, and in- 
terpreted as being due to differences in the possible inter- 
ference pattern - constructive or destructive interference 
- between nonresonant and broad amplitudes. Interfer- 
ence between the nonresonant and /o(1370) amplitudes 
appears to be a plausible explanation for the effect in 
this analysis. As can be seen in Fig. Hj the data exhibit 
a non-trivial interference pattern, and the difficulty in 
modeling this effect leads to model uncertainties in our 
results. 
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TABLE 11: Results of fits to data, with statistical, systematic and model-dependent uncertainties. 
Resonance x y Ax Aj; 

p°(770)7r± 1.0 (fixed) 0.0 (fixed) -0.092±0.036±0.027 io;oi2 0.0 (fixed) 



+0.032 
146 



p°(1450)7r± -0.292±0.071±0.065l2;Jg 0.175±0.078±0.048 0.109±0.080±0.059 IIJ;??^ 0.211±0.073±0.038 

/2(1270)7r± 0.136±0.064±0.040lo;j2g 0.149±0.052±0.030 1":!!?? 0.101±0.063±0.016 io;?g3 -0.248±0.052±0.024 +°:°26 

/o(1370)7r± 0.397±0.067±0.058l!J;;55o -0.151±0.081±0.052 1°;?^^ -0.387±0.064±0.029 IJSl^ga -0.168±0.086±0.055 

Nonresonant -0.200±0.091±0.029 lo;o45 - 0.682±0.070±0.038 to;og2 -0.392±0.089±0.055 io:?28 0.046±0.069±0.055 



TABLE III: Summary of measurements of branching fractions (averaged over charge conjugate states) and CP asymmetries. 
The first error is statistical, the second is systematic and the third represents the model dependence. Also included are 90% 
CL upper limits of the branching fractions of the components that do not have statistically significant fit fractions. 



Mode 


Fit Fraction (%) 


B{B^ Mode)(10-'') 


Acp (%) 


tt^^tt^tt'F Total 




15.2 ± 0.6 ± 1.2 tol 


+3.2 ± 4.4 ± 3.112.0 



p°(770)7r±; p"(770) ^ TT+TT" 53.2 ± 3.7 ± 2.5 1^:^ 8.1 ± 0.7 ± 1.2 l^'.i +18±7±5li. 
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p°(1450)7r±; /9"(1450) TT+TT" 9.1 ± 2.3 ± 2.4 1.4 ± 0.4 ± 0.4 +°-^ 



0.7 

/2(1270)7r±; /2(1270) ^ TT+TT^ 5.9 ± 1.6 ± 0.4 0.9 ± 0.2 ± 0.1 +41±25±13l8^ 

/o(1370)7r±; /o(1370) TT+TT" 18.9 + 3.3 + 2.6^3;^ 2.9 + 0.5 + 0.5 tjjls (< 4.0) +72+15 + 141^ 

7r±7r±7r=F nonresonant 34.9 + 4.2 + 2.9 igj 5.3 + 0.7 + 0.6 1 J'g -lzL+l4 + 7+i7 

/o(980)7r=^; /o(980) ^ Tr+Tr" - < 1.5 

Xcott'*'; XcO tt+tt" - < 0.1 

Xc27r='=; Xc2 TT+TT" - <0.1 



SYSTEMATIC UNCERTAINTIES 

Systematic uncertainties that affect the measurement 
of fit fractions, phases, event yields, and CP asymmetries 
are summarized in Table HVl The fixed SS-background 
yields and asymmetries are allowed to vary and the vari- 
ations of the other fitted parameters are taken as the 
uncertainties. The effect of limited statistics of the data 
sideband and MC samples used to obtain the fixed shapes 
of all the histogram PDFs is accounted for by fiuctuating 
independently the histogram bin contents in accordance 
with their errors and repeating the nominal fit. The un- 
certainties on how well the samples model these distribu- 
tions are also taken into account through various cross- 
checks, including variation of the mass rejection ranges 
and comparison of continuum shapes between sideband 
and signal region in MC samples. 

The fixed parameters of the signal mEs and A_E 
PDFs are studied in the control sample 5+ D'-'tt^; 
D° K^Tr~ . The parameters are determined from data 
and MC samples, from which shift and scale factors are 
calculated and used to adjust the parameters for the nom- 
inal fit. The parameters are then varied in accordance 
with the errors on these shift and scale factors and the 



fits are repeated. Uncertainties due to the mEs distri- 
bution for the qq background, which is fixed in the fit, 
are assessed to be negligible. The fit confirms the value 
of the single parameter of the ARGUS function that is 
taken as input. 

To confirm the fitting procedure, we perform 500 MC 
experiments in which the events are generated from the 
PDFs used in the fit to data. We repeat the exercise 
with qq events alone drawn from the PDF into which 
we embed signal and BB background events randomly 
extracted from the MC samples. Small fit biases are ob- 
served for some of the fit parameters and are included in 
the systematic uncertainties. 

Relative uncertainties in the efficiency, due to PID 
and tracking efficiency corrections are 4.2% and 2.4% re- 
spectively; while Ngj^ has an associated error of 1.1%. 
The efficiency correction due to the selection require- 
ment on the NN output has also been calculated from 
B+ D°TT+; W K+TT- data and MC samples, and 
is found to be (96.2 ± 1.2)%. The error on this correction 
is incorporated into the branching fraction systematic un- 
certainties. 

Measured CP asymmetries could be affected by detec- 
tor charge bias. We include a s yst ematic uncertainty of 
0.005 to account for this effect 21[. Furthermore, some 
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TABLE IV: Absolute maximum values of the systematic uncertainties for the amplitude coefBcients, fit fractions, signal yield, 
and CP asymmetries from various sources described in the text. 



Source 


X 


y 


Aa; 


Ay 


Fit fraction Acp Sif 


;nal yield Signal asymmetry 


BB yields 


0.02 


0.04 


0.02 




0.02 


0.08 


1.4 


0.01 


BB PDF 












0.01 


3.3 




Signal PDF 


0.02 


0.02 


0.01 


0.02 


0.01 


0.08 


48.3 


0.01 


qq Dalitz plot 


0.06 


0.03 


0.05 


0.05 


0.03 


0.14 


47.7 




BB Dalitz plot 


0.03 


0.03 


0.02 


0.02 


0.02 


0.10 


31.6 


0.02 


Efficiency Dalitz 


plot ... 


0.01 








0.03 


0.6 




Fit bias 


0.01 


0.02 


0.02 


0.01 




0.05 


2.8 






- (GeV/cO 



m^.^-(GeV/c') 



FIG. 4: (color online) Dipion invariant mass projection in the 
p°(770) region for (left) B' and (right) _B+ candidates. The 
top row shows all candidates, the middle row shows those with 
cosOh > and the bottom row shows those with cos 9h < 0. 
The colors and shadings follow the same convention as Fig. [3] 



of our selection requirements, for example that on the 
NN output, may induce an asymmetry. We estimate the 
possible size of such an effect as 0.020 based on the study 
of our control sample. 

In addition to the above systematic uncertainties we 
also estimate uncertainties from two sources related to 
the signal Dalitz plot model. The first of these elements 
consists of the parameters of the various components of 
the signal model: the masses and widths of all interme- 
diate resonances, the value of the parameter that char- 
acterizes the nonresonant shape, and the value of the 



Blatt-Weisskopf barrier radius. The associated uncer- 
tainties are evaluated by adjusting the parameters within 
their experimental errors and refitting. The second ele- 
ment is the uncertainty due to the composition of the 
signal model. It reflects observed changes in the param- 
eters of the components when the data are fitted with 
the less significant /o(1370) component removed from the 
model and with one of the states ijj(782), /o(980), XcO, 
or Xc2 added to the model. The uncertainties from each 
of these elements are added in quadrature to obtain the 
final model-dependence. 



DISCUSSION AND SUMMARY 

Our results are shown in Tables [TTl and IIIIl The Dalitz 
plot is dominated by the ^"^(770) resonance and a non- 
resonant contribution which, as seen in other charmless 
three-body hadronic B decays, is well modeled with an 
exponential form-factor. The measured branching frac- 
tion for the decay /9°(770)7r^ agrees with the 
world-average value {sij and is consistent with theoreti- 



cal predictions based on QCD factorization models [44 1 
and SU(3) flavor symmetry The measured branch- 
ing fraction of the nonresonant component is consistent 
with some theoretical predictions [IJ, M,\l^. We find 
the parameter of the nonresonant lineshape to be am- = 
0.28±0.06 GeV~^c* (statistical uncertainties only), which 
is comparable with values found in analyses of other 
charmless decay modes such as B ~> i^TTTr [1^ [l^, [l^ 
and B 



KKK 23, 2 



Contributions from /5°(1450) and /o(1370) are also in- 
cluded in the Dalitz plot model, where the mass and 
width of the /o(1370) are determined to be m^j,(i37o) = 
1400±40 MeV/c2 and r/„(i37o) = 300±80MeV (statistical 
uncertainties only). We have made the first observation 
of the decay B^ /2(1270)7r* with a statistical signifi- 
cance of 6.1f7. After correcting for the /2(1270) — > tt+tt" 
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branching fraction of (84.8_j'2) ^ 10 



^ /2(1270)7r±) 



I 3j| , we obtain 
(24) 



(1.57 ± 0.42 ± 0.16 IJJ:^^) X 10"*' 



where the uncertainties are statistical, systematic, and 
model-dependent, respectively. The latter includes the 
uncertainty of the /2(1270) tt+tt^ branching fraction. 
The above measurements of the branchin g fr actions are 
generally improved from previous results |32l |. although 
some of the uncertainties are not reduced, largely due to 
a more realistic assignment of mo del- dependent uncer- 
tainties in this analysis. 

The 90% confidence-level upper limits for the branch- 
ing fractions of Xcott* and 
to be 



Xc27r* are found 



BiB"^ 



< 1.5 X IQ-^ 

< 2.0 X 10"^ 



(25) 
(26) 



where the Xc(o.2) ^ tt+tt^ widths are determined from 
recent measurements by Belle [i^ . The absence of these 
charmonium contributions precludes the extraction of the 
unitarity triande angle 7 that has been proposed in the 
literature BMSllli. 

We do not find any signal for the decay B^ 
/o(980)7r^. The branching fraction upper limit we obtain 
is consistent with the prediction of a recent perturbative 
QCD calculation if the /o(980) meson is dominated by 



an ss component [47[. 



We do not find any statistically significant CP asym- 
metries for the components in the nominal Dalitz plot 
model. The CP asymmetry in B^ p°(770)7r^ has a de- 
pendence on the presence or absence of the /o(1370) term 
in the model. The CP asymmetry of the /o(1370) term 
itself appears highly sensitive to the Dalitz plot model, 
varying dramatically between the favored and the second 
solution. Since the presence of this component is not es- 
tablished, especially given its insignificant contribution 
in the second solution, we set a 90% CL upper limit on 
its product branching fraction at 4.0 x 10^^. 

In conclusion, we have performed a Dalitz plot analysis 
of B^ TT^TT^TT^ decays based on a 424 fb""'^ data sam- 
ple containing (465 ± 5) x 10^ BB pairs collected with 
the BaBAR detector. Our model includes a momentum- 
dependent nonresonant component and four intermediate 
resonance states: p°(770)7r±, p°{U50)Tr^, /2(1270)7r±, 
and /o(1370)7r±. We do not find any significant contri- 
butions from /o(980)7r='=, Xcott^, or Xc27r^. We find no 
evidence for direct CP violation. Our results will be use- 
ful to reduce model uncertainties in the extraction of the 
CKM angle a from time-dependent Dalitz plot analysis 
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'n^. The results presented here supersede 
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CORRELATIONS BETWEEN FIT FRACTIONS 
AND DIRECT CP ASYMMETRIES 

In Table fVl we present the statistical linear correlations 
between the values of FFj and Acp, j . 

COMPARISON OF RESULTS IN FAVORED AND 
SECOND SOLUTION 

In Table IVTl we give a comparison of the results for the 
two solutions. Note that the Acp value of /o(1370)7r^ is 
at the physical boundary in the second solution. 
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TABLE V: Matrix of statistical correlation coefficients between fit fi'actions and direct CP asymmetries. 





p°(770) 


p''(1450) 


/2(1270) 


/o ( 1 370) Nonresonant 


Parameter 








J^CP 


FF Acp 


FF Acp FF Acp 


p"(770) FF 


l.UU 












p"{770) Acp 


n nn 
u.uu 


l.UU 










p°(1450) FF 


n 1 n 
u.iu 


n 9Q 
u.zy 


l.UU 








p°(1450) Acp 


0.35 


0.14 


0.31 


1.00 






/2(1270) FF 


-0.14 


-0.02 


-0.15 


-0.05 


1.00 




/a (1270) Acp 


-0.12 


0.06 


-0.01 


-0.02 


-0.08 1.00 




/o(1370) FF 


-0.28 


-0.16 


-0.45 


-0.31 


-0.09 0.02 


1.00 


/o(1370) Acp 


-0.14 


0.00 


-0.09 


-0.13 


0.01 -0.10 


-0.14 1.00 


Nonresonant FF 


-0.50 


0.06 


-0.15 


0.06 


-0.17 0.07 


0.26 -0.17 1.00 


Nonresonant Acp 


-0.02 


-0.22 


0.03 


0.09 


0.03 -0.06 


0.13 0.27 0.06 1.00 



TABLE VI: Comparison of results for the favored and second solutions. The uncertainties are statistical only. 



Favored solution 


Second solution 


Inclusive signal yield 


1219 ± 50 


1195 + 45 


Inclusive signal Acp 


+0.032 ± 0.044 


+0.015 + 0.043 


qq background yield 


2337 ± 62 


2358 + 64 


qq background Acp 


+0.002 ± 0.027 


+0.011 + 0.027 




Favored solution 


Second solution 


Resonance 


Fit fraction Acp 


Fit fraction Acp 


p«(770)7r± 


0.532+0.037 +0.18+0.07 


0.458+0.033 +0.03+0.08 


p"(1450)7r± 


0.091+0.023 -0.06+0.28 


0.064+0.016 -0.54+0.24 


/2(1270)7r± 


0.059+0.016 +0.41+0.25 


0.079+0.016 +0.55+0.20 


/o(1370)7r± 


0.189+0.033 +0.72+0.15 


0.030+0.019 -1.00 1^:^^ 


Nonresonant 


0.349+0.042 -0.14+0.14 


0.365+0.042 -0.11+0.14 



